Beyond pressureless gas dynamics: Quadrature-based velocity moment models

Following the seminal work of F. Bouchut on zero pressure gas dynamics which has been extensively used for gas particle-flows, the present contribution investigates quadrature-based velocity moments models for kinetic equations in the framework of the infinite Knudsen number limit, that is, for dilute clouds of small particles where the collision or coalescence probability asymptotically approaches zero. Such models define a hierarchy based on the number of moments and associated quadrature nodes, the first level of which leads to pressureless gas dynamics. We focus in particular on the four moment model where the flux closure is provided by a two-node quadrature in the velocity phase space and provide the right framework for studying both smooth and singular solutions. The link with both the kinetic underlying equation as well as with zero pressure gas dynamics is provided and we define the notion of measure solutions as well as the mathematical structure of the resulting system of four PDEs. We exhibit a family of entropies and entropy fluxes and define the notion of entropic solution. We study the Riemann problem and provide a series of entropic solutions in particular cases. This leads to a rigorous link with the possibility of the system of macroscopic PDEs to allow particle trajectory crossing (PTC) in the framework of smooth solutions. Generalized $\delta$-choc solutions resulting from Riemann problem are also investigated. Finally, using a kinetic scheme proposed in the literature without mathematical background in several areas, we validate such a numerical approach in the framework of both smooth and singular solutions.Comment: Submitted to Communication in Mathematical Science

On the Eulerian Large Eddy Simulation of disperse phase flows: an asymptotic preserving scheme for small Stokes number flows

In the present work, the Eulerian Large Eddy Simulation of dilute disperse phase flows is investigated. By highlighting the main advantages and drawbacks of the available approaches in the literature, a choice is made in terms of modelling: a Fokker-Planck-like filtered kinetic equation proposed by Zaichik et al. 2009 and a Kinetic-Based Moment Method (KBMM) based on a Gaussian closure for the NDF proposed by Vie et al. 2014. The resulting Euler-like system of equations is able to reproduce the dynamics of particles for small to moderate Stokes number flows, given a LES model for the gaseous phase, and is representative of the generic difficulties of such models. Indeed, it encounters strong constraints in terms of numerics in the small Stokes number limit, which can lead to a degeneracy of the accuracy of standard numerical methods. These constraints are: 1/as the resulting sound speed is inversely proportional to the Stokes number, it is highly CFL-constraining, and 2/the system tends to an advection-diffusion limit equation on the number density that has to be properly approximated by the designed scheme used for the whole range of Stokes numbers. Then, the present work proposes a numerical scheme that is able to handle both. Relying on the ideas introduced in a different context by Chalons et al. 2013: a Lagrange-Projection, a relaxation formulation and a HLLC scheme with source terms, we extend the approach to a singular flux as well as properly handle the energy equation. The final scheme is proven to be Asymptotic-Preserving on 1D cases comparing to either converged or analytical solutions and can easily be extended to multidimensional configurations, thus setting the path for realistic applications

Sensor validation of a Structural Health Monitoring Process for Aircraft Nacelle

This paper details the implementation process of an embedded structural health monitoring (SHM) system enabling condition-based maintenance of aircraft nacelles. One critical issue before being able to make use of such system is to ensure the effective bonding of the chosen actuators and sensors with their host structure, especially as the latter will be exposed to harsh environments and wide operational variability. In this work, we are concerned with the composite components of the nacelle and we use piezoelectric elements as both sensors and actuators. We propose an integrated approach that allows to validate a combination “Substrate—Glue—Piezoelectric” (SGP) and thus provides criteria to choose and size these assemblies. This validation scheme is based on the observation of the variations of the static capacity of the piezoelectric element after enduring various temperature and stress conditions when bonded to its host structure. Based on those SGP combinations, an active SHM strategy interrogating the structure by means of elastic wave propagation is currently being developed and preliminary results on samples representative of the nacelle are presented and discussed.Projet AIRCELLE (EPICE/CORALIE

Thermoacoustic instability - a dynamical system and time domain analysis

This study focuses on the Rijke tube problem, which includes features relevant to the modeling of thermoacoustic coupling in reactive flows: a compact acoustic source, an empirical model for the heat source, and nonlinearities. This thermo-acoustic system features a complex dynamical behavior. In order to synthesize accurate time-series, we tackle this problem from a numerical point-of-view, and start by proposing a dedicated solver designed for dealing with the underlying stiffness, in particular, the retarded time and the discontinuity at the location of the heat source. Stability analysis is performed on the limit of low-amplitude disturbances by means of the projection method proposed by Jarlebring (2008), which alleviates the linearization with respect to the retarded time. The results are then compared to the analytical solution of the undamped system, and to Galerkin projection methods commonly used in this setting. This analysis provides insight into the consequences of the various assumptions and simplifications that justify the use of Galerkin expansions based on the eigenmodes of the unheated resonator. We illustrate that due to the presence of a discontinuity in the spatial domain, the eigenmodes in the heated case, predicted by using Galerkin expansion, show spurious oscillations resulting from the Gibbs phenomenon. By comparing the modes of the linear to that of the nonlinear regime, we are able to illustrate the mean-flow modulation and frequency switching. Finally, time-series in the fully nonlinear regime, where a limit cycle is established, are analyzed and dominant modes are extracted. The analysis of the saturated limit cycles shows the presence of higher frequency modes, which are linearly stable but become significant through nonlinear growth of the signal. This bimodal effect is not captured when the coupling between different frequencies is not accounted for.Comment: Submitted to Journal of Fluid Mechanic

Task-based adaptive multiresolution for time-space multi-scale reaction-diffusion systems on multi-core architectures

A new solver featuring time-space adaptation and error control has been recently introduced to tackle the numerical solution of stiff reaction-diffusion systems. Based on operator splitting, finite volume adaptive multiresolution and high order time integrators with specific stability properties for each operator, this strategy yields high computational efficiency for large multidimensional computations on standard architectures such as powerful workstations. However, the data structure of the original implementation, based on trees of pointers, provides limited opportunities for efficiency enhancements, while posing serious challenges in terms of parallel programming and load balancing. The present contribution proposes a new implementation of the whole set of numerical methods including Radau5 and ROCK4, relying on a fully different data structure together with the use of a specific library, TBB, for shared-memory, task-based parallelism with work-stealing. The performance of our implementation is assessed in a series of test-cases of increasing difficulty in two and three dimensions on multi-core and many-core architectures, demonstrating high scalability

High order time integration and mesh adaptation with error control for incompressible Navier-Stokes and scalar transport resolution on dual grids

International audienceRelying on a building block developed by the authors in order to resolve the incompressible Navier-Stokes equation with high order implicit time stepping and dynamic mesh adaptation based on multiresolution analysis with collocated variables, the present contribution investigates the ability to extend such a strategy for scalar transport at relatively large Schmidt numbers using a finer level of refinement compared to the resolution of the hydrody-namic variables, while preserving space adaptation with error control. This building block is a key part of a strategy to construct a low-Mach number code based on a splitting strategy for combustion applications, where several spatial scales are into play. The computational efficiency and accuracy of the proposed strategy is assessed on a well-chosen three-vortex simulation

On the Direct Numerical Simulation of moderate-Stokes-number turbulent particulate flows using Algebraic-Closure-Based and Kinetic-Based Moment Methods

In turbulent particulate flows, the occurrence of particle trajectory crossings (PTC) is the main constraint on classical monokinetic Eulerian methods. To handle such PTC, accounting for high-order moments of the particle velocity distribution is mandatory. In the simplest case, second-order moments are needed. To retrieve these moments, two solutions are proposed in the literature: the Algebraic-Closure-Based Moment Method (ACBMM) and the Kinetic-Based Moment Method (KBMM). The ACBMM provides constitutive relations for the random-uncorrelated-motion (RUM) particle kinetic stress tensor as algebraic closures based on physical arguments (Simonin et al. 2002; Kaufmann et al. 2008; Masi 2010; Masi & Simonin 2012). These closures rely on the internal energy, namely the RUM particle kinetic energy, which is obtained using an additional transport equation. Alternatively, it is possible to directly solve for the second-order moment by providing a closure for the third-order correlation. The KBMM proposes a kinetic description, that is, the number density function (NDF) is reconstructed based on the resolved moments and on a presumed shape. In the present work, an isotropic Gaussian and the anisotropic Gaussian closure of Vié et al. (2012) are used. The goal of the present study is to provide a first comparison between ACBMM and KBMM, using the same robust numerical methods, in order to highlight differences and common points. The test case is a 2D turbulent flow with a mean shear

Turbulent thermal convection driven by heated inertial particles

The heating of particles in a dilute suspension, for instance by radiation, chemical reactions or radioactivity, leads to local temperature fluctuations in the fluid due to the non-uniformity of the disperse phase. In presence of a gravity field, the fluid is set in motion by the resulting buoyancy forces. When the particle density is different than the fluid, the fluid motion alters the spatial distribution of particles and possibly strengthens their concentration inhomogeneities. This in turn causes more intense local heating. Direct numerical simulations in the Boussinesq limit show this feedback loop. Various regimes are identified depending on the particle inertia. For very small particle inertia, the macroscopic behavior of the system is the result of many thermal plumes that are generated independently of each other. For significant particle inertia, clusters of particles are observed and their dynamics control the flow. The emergence of very intermittent turbulent fluctuations shows that the flow is influenced by the larger structures (turbulent convection) as well as by the small-scale dynamics that affect particle segregation and thus the flow forcing. Assuming thermal equilibrium between the particles and the fluid (i.e. infinitely fast thermal relaxation of the particle), we investigate the evolution of statistical observables with the change of main control parameters (namely the particle number density, the particle inertia and the domain size), and propose scaling argument for these trends. Concerning the energy density in the spectral space, it is observed that the turbulent energy and temperature spectra follow a power law, the exponent of which varies continuously with the Stokes number. Furthermore the study of the spectra of the temperature and momentum forcing (and thus of the concentration/temperature and velocity/temperature correlations) gives strong support to the proposed feedback loop mechanism. We then discuss the intermittency of the flow, and analyze the effect of relaxing some of the simplifying assumptions, thus assessing the relevance of the original studied configuration